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Abstract 

The Black-Scholes formula for pricing options on stocks and other securities has been gener- 
ahzed by Merton and Garman to the case when stock volatihty is stochastic. The derivation 
of the price of a security derivative with stochastic volatihty is reviewed starting from the 
first principles of finance. The equation of Merton and Garman is then recast using the path 
integration technique of theoretical physics. The price of the stock option is shown to be the 
analogue of the Schrodinger wavefuction of quantum mechacnics and the exact Hamiltonian 
and Lagrangian of the system is obtained. The results of Hull and White are generalized to 
the case when stock price and volatility have non-zero correlation. Some exact results for 
pricing stock options for the general correlated case are derived. 
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I Introduction 



The problem of pricing the European call option has been well-studied starting from the 
pioneering work of Black and Scholes [1] , [2] . The results of Black and Scholes were general- 
ized by Merton, Garman [3] and others for the case of stochastic volatility and they derived 
a partial differential equation that the option price must satisfy. 

The methods of theoretical physics have been applied with some success to the problem 
of option pricing by Bouchaud et al [4]. Analyzing the problem of option pricing from the 
point of view of physics brings a whole collection of new concepts to the field of mathematical 
finance as well as adds to it a set of powerful computational techniques. This paper is a 
continuation of applying the methodology of physics, in particular, that of path integral 
quantum mechanics, to the study of derivatives. 

It should be noted that, unhke the paper by Bouchaud and Sornette [4] in which the 
pricing of derivatives is obtained by techniques which go beyond the conventional approach 
in finance, the present paper is based on the usual principles of finance. In particular, a 
continuous time random walk is assumed for the security and a risk- free portfolio is used 
to derive the derivative pricing equation; these are reviewed in Sec II. The main focus of 
this paper is to apply the computational tools of physics to the field of finance; the more 
challenging task of radically changing the conceptual framework of finance using concepts 
from physics is not attempted. 

An explicit analytical solution for the pricing of the European call (or put) option for 
the case of stochastic volatility has so far remained an unsolved problem. Hull and White 
[5] studied this problem and obtained a series solution for the case when the correlation of 
stock price and volatility, namely r, is equal to zero, i.e. p = (uncorrelatcd). Hull and 
White [5] also obtained an algorithm to numerically evaluate the option price using Monte 
Carlo methods for the case of p 7^ 0. 

Extensive numerical studies of option pricing for stochastic volatility based on the algo- 
rithm of Finucane [6] have been carried out by Mills ct al [7] . 

The main focus and content of this paper is to study the problem of stochastic volatility 
from the point of view of the Feynman path integral [9] . The Feynman-Kac formula for the 
European call option is well known [2, 8]. There have also been some applications of path 
integrals in the study of option pricing [10, 11]. In this paper the path integral approach to 
option pricing is first discussed in its generality and then applied to the problem of stochastic 
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volatility. 



The advantage of recasting the option pricing problem as a Feynman path integral is 
that this allows for a new point of view and which leads to new ways of obtaining solutions 
which are exact, approximate as well as numerical, to the the pricing of options. 

In Section II the differential equation of Merton and Garman [3] is derived from the 
principles of finance. In Section III this equation is recast in the formalism of quantum 
mechanics. In Section IV a discrete time path integral expression is derived for the option 
price with stochastic volatility which generalizes the Feynman-Kac formula. In Section V 
the path integration over the stochastic stock price is performed explicitly and a continuous 
time path integral is then obtained. And lastly in Section VI some conclusions are drawn. 



We review the principles of finance which underpin the theory of security derivatives, and in 
particular that of the pricing of options. We will derive the results not using the usual method 
used in theoretical finance based on Ito-calculus (which we will review for completeness), 
but instead from the Langevin stochastic differential equation. Hence we start from first 
principles. 

A security is any financial instrument which is traded in the capitals market; this could 
be the stock of a company, the index of a stock market, government bonds etc. A security 
derivative is a financial instrument which is derived from an underlying security and which 
is also traded in the capitals market. The three most widely used derivatives are options, 
futures and forwards; more complicated derivatives can be constructed out of these more 
basic derivartives. 

In this paper, we analyze the option of an underlying security S. A European call option 
on S is a financial instrument which gives the owner of the option to buy or not to buy the 
security at some future time T ^ t for the strike price of K. 

At time t — T, when the option matures the value of the call option f(T) is clearly given 



II Security Derivative with Stochastic Volatihty 



by 




S{T) > K 
S{T) < K 



(la) 



(16) 
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The problem of option pricing is the following: given the price of the security S{t) at 
time t, what should be the price of the option / at time t < T7 Clearly / = f{t, S{t), K, T). 
This is a final value problem since the final value of / at i = T has been specified and its 
initial value at time t needs to be evaluated. 

The price of / will be determined by how the security S (t) evolves in time. In theoretical 
finance, it is common to model S{t) as a random (stochastic) variable with its evolution 
given by the stochatic Langevin equation (also called an Ito- Wiener process) as 

^ (ps{t) + aSR{t) (2a) 

(JjL 

where R is the usual Gaussian white noise with zero mean; since white noise is assumed to 
be independent for each time t, we have the Dirac delta-function correlator given by 

< RtRf >= 5{t-t'). {2b) 

Note is the expected return on the security S and a is its volatihty. White noise R{t) 
has the following important property . If we discretize time t — ne, then the probability 
distribution function of white noise is given by 



P(Rt) = V^^"' ' (3a) 
For random variable R it can be shown using eqn. (3a) that 

Rt — — h random terms of 0(1) (36) 

In other words, to leading order in e, the random variable i?| becomes deterministic. This 
property of white noise leads to a number of important results, and goes under the name of 
Ito calculus in probability theory. 

As a warm-up, consider the case when a — constant; this is the famous case considered 
by Black and Scholes. We have 

^ ^ li^ f{t + e,S{t + e))-f{t,S{t)) 

or, using Taylors expansion 

d^^dl Old_S ld^d_S 

dt dt dS dt 2dS^^ dt^ ^ ' ^ ' 
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But 

fi2 = a^S^R^ + 0(1) 

CLv 



= -^a^S^ + 0(1) (4c) 
Hence we have, for e — > 

Since eqn. (4d) is of central importance for the theory of security derivatives we also give 
a derivation of it based on Ito-calculus. Rewrite eqn. (2a) in terms differentials as 

dS = 4>Sdt + aSdz 

where the Wiener process dz — Rdt, with R{t) being the Gaussian white noise. Hence from 
eqn. (3b) 

{dzf = Rlidtf ^dt + 0{dt^/'') 

and hence 

{dSf = a^S^dt + 0{df/^). 
Prom the equations for dS and {dS^ given above we have 

and using ^ = i? we obtain eqn. (4d). 

The fundamental idea of Black and Scholes is to form a portfolio such that, instan- 
taneously, the change of the portfolio is independent of the white noise R. Consider the 
portfolio 

df 

^ = f- Q^S (5a) 

i.e. TT is a portfolio in which an investor sells an option / and buys || amount of security 
5". We then have from eqns. (4b) and (5a) 
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dn df df dS 
dt dt dS dt 

df 1 d'f 

(56) 

The portfolio tt has been adjusted so that its change is deterministic and hence is free 
from risk which comes from the stochastic nature of a security. This technique of canceUing 
the random fluctuations of one security (in this case of /) by another security (in this case 
S) is called hedging. The rate of return on hence must equal the risk-free return given by 
the short-term risk- free interest rate r since otherwise one could arbitrage. Hence 

— = m (oc) 
dt ^ ' 

which yields from eqn.(5b) the famous Black-Scholes [1] equation 

Note the parameter of eqn. (2a) has dropped out of eqn. (5d) showing that a risk- 
neutral portfolio TT is independent of the investors expectation as reflected in 0; or equiv- 
alently, the pricing of the security derivative is based on a risk-neutral process which is 
independent of the investors opinion. 

We now address the more complex case where the volatility a{t) is considered to be 
a stochastic (random) variable. We then have the following coupled Langevin equations 
(following Hull and White [5]), for cr^ = V, 

— = (j)S + aSR (6a) 

(Jjv 

dV 

— = ^V + ^VQ (66) 

where n is the expected rate of increase of the variance V, ^ its volatility and R h Q are 
correlated Gaussian white noise with zero means and with correlators 

< RtRf >=< QtQt' >= S{t - t') (6c) 

and 

< RtQt' >= p5{t - t') {6d) 
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with < 1 being the reflection of the correlation between S and V. 

Note that one can choose a different process from eqn. (6) for stochastic volatility as has 
been done by various authors [5] but the main result of the derivation is not affected. 

The procedure for deriving the equation for the price of an option is more involved when 
volatility is stochastic. Since volatility V is not a traded security this means that a portfolio 
can only have two instruments, namely / and S. To hedge away the all fluctuations due to 
both S and V, we need at least three flnancial instruments to make a portfolio free from the 
investors subjective preferences. Since we have only two instruments from which to form our 
portfolio we will not be able to completely avoid the subjective expectations of the investors. 

Consider the change in the pricing of the option in the presence of stochastic volatility; 
using eqns (6a) and (6b) we have from Taylors expansion 



% = ]imhf{t + e,V{t + e),S{t + e)-f{t,V{t),Sm} 

at e^o e 



+ 



Using 



1 / ^l^^2 ^ ^l^^2 ^ o d\f dVdS 



2\dS^Ut^ dV^^dt^ dSdV dt dt 



Rt = Qt = +0(1), 
RtQt =- + 0(1) 



we have, similar to eqn. (4a), for e — > 

df Of ^^df ^df 1 
Jt = Tt^'^^dS^^'^dV^2 



a2f a2f 



dVdS 



Define the deterministic and stochastic terms by 

df 

^ = + faiR + fa2Q 



(7a) 



(76) 



(7c) 



{7d) 



Consider a portfoho which consists of 9i{t) amount of option / and 92{t) amount of stock 
S, that is 

7r{t)^0i{t)f{t) + 02{t)S{t) (8a) 
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Then 



— = 01— + 02— + —f + —S 
dt dt dt dt dt 



m 



We consider a self-replicating portfolio in which no cash flows in or out and the quantities 
Oiit) and 92{t) are adjusted accordingly. It then follows that 



Hence 



Itt 



dn „ df „ dS 
— = 9i— + 92 — 
dt dt dt 



9ifu + (t)92S + (aii? + a2Q)9if + 92SaR 



We choose the portfolio such that, in matrix notation 

CKi (T 9if 

= 

q;2 J [ 92S 

which yields 



d-K 
'dt 



9ifuj + 4)92S 



(8c) 



(8e) 



(9a) 



Note that all the random terms from the rhs of eqn. (8d) have been removed by our choice 
of portfoho which satisfies the constraint given in eqn. (8e). Since ^ is now deterministic, 
the absence of arbitrage requires 

{%) 



dn 
'dt 



or, from eqns. (8d) and (9a) 



{u; - r)9J + {(f)- r)92S = 



(9c) 



We hence have from eqns. (8e) and (9c) that eqn. (9c) above can be satisfied in general by 
making the following choice uj and 0, namely 



0L2 



[a;-r,0-r] = [Ai(t),A2(t)] 
where Ai and A2 are arbitrary. In components, we have from above 

LO — r — Aicci + A2Q;2 
4> — r = Aicr 



m 

(9e) 
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Eliminating Ai yields using eqn. (9e) yields from eqn. (9d) 



fu-rf ^ 



+ A2/q;2 



and hence from eqns. (7c) and (7d) we obtain the Merton and Garman [3] equation given 

by 



dt 2 



dSdV 



-rf = -rS^ + [X2 - 1^] (10a) 



Note the fact that both and Ai have been eliminated from (10a) since we could hedge 
for the fluctuations of S by including it in the portfolio vr. In order words the subjective 
view of the investor visa vis security S has been removed. However, the appearance of A2 in 
eqn. (10a) reflects the view of the investor of what he expects for the volatility of S. Since 
V is not traded, the investor can't execute a perfect hedge against it as was the case for 
constant volatility and hence the pricing of the option given in eqn. (10a) has a subjective 
factor in it. 

It has been argued by Hull and White [5] that for a large class of problems, we can 
consider 1 to be a constant which in effect simply redefines // to be — A2 . 
Hence we have 



dt 2 



2a2*^J I So -' 



9v 



rf--rS%-,a^% (106) 



Eqn. (10) above is an equation for any security derivative with stochastic volatility. 
Whether it is an European option or an Asian option or an American option is determined 
by the boundary conditions that are imposed on /; in particular the boundary condition 
given in eqn. (1) is that of an European call option. 



Ill Black-Scholes-Schrodinger Equation for Option Price 

We recast the Merton and Garman equation given in eqn. (10) in terms of the formalism 
of quantum mechanics; in particular we derive the Hamiltonian for eqn. (10) which is the 
generator of infinitesimal translations in time. Since both S and V are random variables 
which can never take negative values, we define new variables x and y 
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S = e"", -oo < a; < oo (11a) 
V = e^,-oo<y<oo. (116) 

Prom cqn. (10b), we then obtain the Black-Scholes-Schrodinger equation for the price of 

security derivative as 

% = (H + r)f (12a) 

where / is the Schrodinger wave function. The "Hamiltonian' H is a differential operator 
given by 

Prom eq. (12a) we have the formal solution 

/ = e*^^^'-)/ (13) 

where / is fixed by the boundary condition. Note from eqn. (13) that option price / seems 
to be unstable being represented by a growing exponential; however, as will be seen later 
since the boundary condition for / is given at final time T eqn. (13) will be converted to a 
decaying exponential in terms of remaining time = T — t. 

Let p{x, y, r\x', y'), be the conditional probability that, given security price x' and volatil- 
ity y' at time r = 0, it will have a value of x and volatility y at time t. We have the boundary 
condition at r = given by Dirac delta-functions, namely 

p{x, y, 0\x', y') = 5{x' - x)5{y' - y) (14) 

The derivative price is then given, for < i < T, by the Peynman-Kac formula Refs [5, 
8] as 

f{t; X, y) = e-(^-*)'- / dx'p{x, y,T- tx')g{x') (15a) 

J —CO 

where 

/ + 00 
dy'p{x,y,T -tx',y') (156) 
-oo 

Note that integration over y' is decoupled from the function g{x') due to eqn. (1). It follows 
from eqn. (14) that f{t;x,y) given in (15a) satisfies the boundary condition given in eqn. 
(1), i.e. 
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f{T,x,y)^g{x) (156) 

We rewrite eqns. (15a) and (15b) in the notation of quantum mechanics. The function 
f{x) can be thought of as an infinite dimensional vector \f > of a function space with 

components f{x) =< x\f > and f*{x) =< f\x> (where * stands for complex conjugation). 
For — oo < X < oo the bra vector < a; > is the basis of the function space and the ket vector 
|a; > is its dual with normalization 

<x\x' >^5{x-x'). (16) 
The scalar product of two functions is given by 
<f\g> = dxf*{x)g{x) 

/OO r 
dx\x >< x\}\g > 

and yields the completeness equation 

/oo 
dx\x >< x\ (17a) 
-oo 

where / is the identity operator on function space. For the case of stochastic stock price and 
stochastic volatility, we have 

/oo 
dxdy\x,y >< x,y\ (176) 
-oo 

where \x,y >= \x > <^\y >. 

From eqn (13) we have, in Dirac's notation 

|/,i>=e*(^+^)|/,0> 

and boundary condition given in eqn. (lb) yields 

|/,T>=e*(^+^)|/,0>= \g>. 

Hence 

>=e(^-*)(^+^)|^> (18a) 

or, more explicitly 

f{t,x,y) = <x,y\f,t> 

= e-'-M <x,y|e-M^b> (186) 
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or, using completeness equation (17b) 

fit, X, y) = e-'^^-*) r dx/dy' < x, y\e-^^-^^^ \x' , y' > g{x') (19a) 

J—oo 

which yields for remaining time t — T — t, from eqn. (15b) 

p{x, y, t\x', y') =< X, y|e-^-^|a;', y' > . (196) 

Note remaining time runs backwards, i.e. when r = 0, we have r = T and when t = T, real 
time t — 0] as mentioned earlier, expressed in terms of remaining time r the derivative price 
is given by a decaying exponential as given in eqn. (19b). 

Before tackling the more complicated case of option pricing with stochastic volatility I 
apply the formalism developed so far to the simpler case of constant volatility; the detailed 
derivation is given in Appendix 1. The well known results of Black and Scholes are seen to 
emerge quite naturally in this formalism. 



IV. The Discrete-Time Feynman Path Integral 

The Feynman path integral [9] is a formulation of quantum mechanics which is based on 
functional integration. In particular, the Feynman path integral provides a functional inte- 
gral realization of the conditional probability p{x, y, t\x', y'). To obtain the path integral we 
discretize time r into N points with spacing e = t/N. Then, from (19b), 

p{x, y, t\x', y') = lim < x, yje"'^ • • • e-'"\x', y' > . (20) 

A"— »oo 

Hence forth, we will always assume the — > oo limit. 

Inserting the completeness equation for x and y, namely eqn (17b), (A^ — 1) times in 
eqn. (20) yields 

iV-l jV 

p{x, y, r\x', y') = {'[[ / dxidyi) [| < Xi, yi|e~'^|xi_i, y^-i > (21a) 
with boundary conditions 

xq = x\ yo = y' (216) 
xn ^x, yo = y- (21c) 
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We show in the Appendix 2 that for the H given by eqn (12b) we have the Feynman 
relation given in (A9) and (A13) 

< Xi, yi\e-'^\xi^i, >= ——==e'^^'> (22a) 

27re ^ VI - P 

where the 'lagrangian' L is given in eqn (A15), for 6x = Xi — Xi-i and dy = yi — yi-i as 

For e — > 0, we have, as expected 

< Xi, yi\e~'"\xi-i, yi-i >= d{xi - Xi-i)d{yi - yi-i) + 0(e) (23) 

and the prefactors to the exponential on the rhs of (22a) ensure the correct limit. 
We form the 'action' S by 

N 

S = eJ2m + 0ie). (24) 



i=l 

Note S is quadratic in Xi and non-linear in the y^ variables. 

Hence, from eqns. (20), (21), (22) and (24) we have the path integral 

/oo 
dy'p{x,y,T\x',y') (25a) 
-oo 

= lim / DXDYe^ (256) 

where, lor e — t/N 

DX^ ' U (25c) 

-I ^27re(l - p2) t\ J-^ ^27re(l - p^) 

and 



•t=i 



Note we have included an extra dy^ — dy' integration in DY given in (25d) due to the 
integration over y' in eqn (25a). 

Eqn (25) is the path integral for stochastic stock price with stochastic volatility in its 
full generality. To evaluate expressions such as the correlation of S{t) with V{t'), the path 
integral in eqn (25) has to be used. 
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Since the action S is quadratic in the Xi, the path integral over the Xi variables in eqn. 
(25) can be done exactly. The details are given in Appendix 3. 
Prom eqn (A21) we obtain 

/gS'o+S'i 
BY (26a) 
V2^re(1^7)EE^ 

where the first term in eqn. (22b) gives 

S.= -^i:(f + .4)^ (266) 

and 5"! is the result of the DX path integration given by eqn (A21b) as 

1 ^ 1 

2(1 -p2)eEi=ieS'' ^ 2 

-^^e^^/2[% + 6(/.-yF. (26c) 

Extensive numerical studies of the pricing of European call option has been carried out 
in [13] based on eqns. (26a-c). 



V. Continuous Time Feynman Path Integral 

We first derive the continuum hmit of the Black-Scholes constant volatihty case before an- 
alyzing the more complex case of stochatic volatihty. The DXbs - path integral for the 
Black-Scholes case , from eqn (A5b), has a measure DXbs = [ 27rL^ ]^^^ HiLi which is 
essentially the measure for the flat space R^; we can hence take the — * oo for DX and 
obtain a well defined continuous-time path integral [9]. From eqn (A5a) we obtain taking 
the continuum limit of e — > 

with boundary conditions x(0) = x' and x{t) — x. We see from above that the Black-Scholes 
case of constant volatility corresponds to the evolution of a free quantum mechanical particle 
with mass given by l/cr^ . This yields the Black-Scholes result 

Pbs{x,t\x') =< x|e-^^^^|x' [ DXBse^"^ 
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The exact result for Pbs{x,t\x') is given in eqn (A4). 

Eqn. (26) provides a mathematically rigorous basis for taking the N ^ oo hmit for the 
case of stochastic volatility. On exactly performing the / DX-path integral, the remaining 
/ DY - path integral has a measure DY — [j^^]^^'^ HfLi dVi which, just as in the Black- 
Scholes case, is essentially the measure for the flat space (unhke the nonlinear expression 
in eqn. (25c)); we can hence take the — > oo for DY and for 5*0 + Si and obtain a well 
defined continuous-time path integral. 

We take the hmit of e — > 0; we have, t — ie, eJ^iLi e^* — ^ lo dte^^*^ = tw, and \5yi — > ^. 
Hence, from (26) we have 

S= So + Si (27a) 



2(1 -p^ W*- 2 Jo 



2{1 -p'^)w 



with boundary value 



and 



li^-\)£dten^yr (27c) 
y{r) = y {27d) 



p{x, y, t\x') = f DY / . (28) 

■' V 27r(l — p^)TW 



Note taking the continuum hmit is possible only after the discrete / DX path integration 
has been performed since the nonlinear measure given by (25c) does not have a finite con- 
tinuum limit. All correlation functions < S{t)V{t') > can be obtained from the continuum 
limit of the discrete correlators < e^^e^™ >. 

For the case of p = 0, to recover the Black-Scholes formula, set w = a"^, where a = e^^^ 
is the volatility at time t = 0. As has been noted by Hull and White [5], for the case of 
p = the conditional probability depends on stochastic volatility y{t) only through the 
combination w — ^ Jq e^^'^^/^dt. For p 7^ 0, we see from eqn. (27c) that p{x, y, tx') now 

depends on w as well as on m = - / e^^^^^'^dt and e^^^^/^. We hence have from (27) and 

r Jo 

(28) 

POO qSi{u,v,w) 

p{x,y,T\x') = / dw du dv . g{u,v,w) (29a) 

•^0 ^27r(l - p2)™ 
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where 

Si{u, V, w) = -— — ^ {x -x' + rT-^w+ '^{v - e^/^) 

2(1 — p^)tw 2 ^ 

-^(/^-yM^- (296) 

and g{u, v, w) is given by the path integral 

g(u,v,w) = / DYe^°S{w - - f e''^'^dt}5{v - e^(^)/'} 

J T Jo 

5{u-- re^^'^/^dt} (29c) 

T Jo 

Prom eqn (29c) we see that g{u, v, w) is the probabihty density for u, v, and w. The path 
integral for g{u, v, w) is nonlinear and cannot be performed exactly. We see that p{x, y, t\x') 



is the weighted average of the integrand e'^^ Vv27r(l — p^)tw with respect to g{u, v, w). 

Note g{u,v,w) has the remarkable property that it is independent of p. Following Hull 
and White [5] one can expand the integrand in (29a) in an infinite power series in u,v,w 
and reduce the evaluation of p{x,y,T\x') to finding all the moments of u,v and w; in other 
words we need to evaluate 

<u''w"'vP>= / du dw v!" w"^ g{u,v,w) (30a) 

= / DY[- /%^(*)/2d^]nrl f\Yit)^nm^YiO)/2^So (gg^,) 
J T Jo T Jo 

The path integral for < u^w^^v^ > can be performed exactly. Rewrite (30b) as 

< U'^W'^VP >= ^ dti--- dtndtn+l ■ ■ ■ dtn+mZ{j, V, p) (30c) 

where 

Z{j,y,p) = J DYeM[dtj{t)y{t)]eP'^^'y'e'° (31a) 
and from eqns. (30b) and (30c) we have 

m = ^ E ^(^ - ^0 + E ^(^ - 1^) + Im (316) 

^ i=l i=n+l ^ 

n+m 

= E ai5(t-U) + ^5(t). (31c) 
The path integral for Z{j,y,p) is evaluated exactly in Appendix 3 and yields 

< uW >= a^^^ n / dU^' (32a) 
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where from eqns (A41a), (A41b) and (31c) we have, after some simphfications 

-j^ n+m n+m 

F^if^- -e') E ("iti + E - tj)^! + F' (326) 

with 

F' = ^Pr(/. - le^) + \pe E^ + \p'er. (32c) 

Note that in obtaining eqn (32) wc have used the identity that ^^(0) = |, as given in 
eqn. (A37a) and reflects the identity Jq dtd(t — r) = |. AU the tj integrations in (32) 
can be performed exactly since the exponent is hnear in the tj 's. The expression (32) for 
< u^w^v^ > generahzes the result of Hull and White [5] as we have an exact expression for 
all the moments of u, v and w as well as for their cross-correlators. 

We explicitly evaluate the first few moments using eqn (32). 

< w >= — / dte^* 



T jQ 

— {e^^ - 1) (33) 

^1T 



where V — 



2V- ' 



r2 



+ 



.(e+^)(2/x+a M/i+a M2/^+a 

We evaluate < w'^ > for the case of = 0; we have from eqn (32) 



(34) 



< >-- 



T'* Jo Jo Jo 



(e^^'^ - 9e«'^ + 8 + eeV) (35) 



3^6x2 

Eqns (33), (34) and (35) agree exactly with the results stated (without derivation) in 
Hull & White [5] . It is reassuring to see that two very different formalisms agree exactly and 
this increases ones confidence in the path integral approach. 

We compute a few moments that have not been given in Hull and White [5] . We further 
have (recall a — e^/^ ) 
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Y 11 

T Jo 



2a 



g(M/2-^V8)r _ 



(36a) 
(366) 



,(p/2-e/8)T 



2eW2-€V4)r 



_ (/.2 + e/4) (/X - e/8) (/X + e/4) (/. - em -em){^- e/s) 

and lastly, for oi = |, 02 = 1 in eqn (31c), we have 



< uw >- 



T"* 70 70 



g(M/2-^V8)r 



+ 



(37) 



[3(/. + e/2)(/x + e/4) (/^ + e/2)(/.2 - e/4) ' 3(/. - evi6) J ■ ^^^^ 

All the other moments < u'^w^^v^ > can similarly be evaluated, which in turn yields an 
infinite series solution for p{x,y,T\x') in eqn (29a). 



VI. Conclusion 

The complete information regarding the dynamics of how stock price S{t) and its volatility 
V{t) evolve, their cross-correlators as well as their fluctuations is given by the discrete time 
path integral given in eqn. (25). In particular, to numerically study correlators such as 



> wc need to start from the path integral given in eqn. (25). 



However, if one is interested solely in the price of the option, one needs to only determine 
p{x, y, t\x', y') and in this case the simplified discrete-time path integral obtained in eqn. (26) 
should be used as the starting point for numerical studies. A detailed numerical simulation 
based on the results of this paper has been carried out in rcf [13]; in particular, it has been 
shown that implied volatility smile and frown can be obtained for different values of the 
correlation parameter p as well as for varying values of S{t). 

In summary, we reformulated the option pricing problem in the language of quantum 
mechanics. We then obtained a number of exact results for the option price of a security 
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derivative using the path integral method. Stochastic volatihty introduces a high order of 
nonhnearity in the option pricing problem and needs to be studied using approximate and 
numerical techniques. 

The continuous time path integral given in eqn. (29c) is a nonlinear path integral which 
has no known exact solution. Various techniques could be used to obtain new approximations 
for g(u,v,w) relevant for special applications. 

Numerical algorithms based on the path integral provide a wide range of new numerical 
algorithms which go beyond the usual binomial tree and its generalizations. 

The comparative efficiency of the different algorithms is being studied. The Feynman 
path integral has been numerically studied using many algorithms and which could prove 
useful in addressing problems in finance. 
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Appendix 1 Black-Scholes Formula for Constant Volatil- 
ity 

For the case of constant volatility the Black-Scholes equation is given by the following [2] 

dt~ 2" ^ ds^ ^^dS^"^^ ^^^^ 
Making the change of variable as in (2) 

S = e^, —00 < X < 00 

yields 

with the Black-Scholes Hamiltonia given by 

To evaluate the price of the European call option with constant volatility we have the 
Feynman-Kac formula 

/CO 
dx'<x|e-M^««|x'>^(x') 
-00 

To compute Pbs{-i^,tW) =< xlc^^^'^^'lx' >, where remaining time is r = T — t, we use 
the Hamiltonian approach to illustrate another powerful technique of quantum mechanics. 
We go to the 'momentum' basis in which Hbs is diagonal. The Fourier transform of the 
\x > basis to 'momentum space' is given by 

< x\x' >= 6{x - x') - /" ^e*f(^-^') 

dp 

— / — < x\p >< p\x > 
J -00 2tt 

which yields for momentum space basis \p > the completeness equation 

i^r^\p><p\ (A3) 

and the scalar product 



< x\p >= e*^""; < p\x >= e" 



tpx 
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Prom (A2) and the equation above we have the matrix elements of H given by 

< x\Hbs\p >= HBse'"'' = + ^(^^' - r)p}e^''^ 

Using (A3) and the equation above yields 
< x\e-^^''^\x' > = f°° ^ x|e"^-^^«|p >< p\x' > 

J-oo 2% 



— oo 27r 2 

■■exp[--^{x-x' + T{r-ay2)y] (A4) 



' 2t(72 

The result above is the Black-Scholes distribution. Recall x' — log{S{T)),x — log(5'(i)) 
and T — T — t; the stock price evolves randomly from its given value of S{t) at time t to 
a whole range of values for S{T) at time T. Equation (A4) above states that log(5'(T')) 
has a normal distribution with mean equal to log{S{t)) + (r — (7^/2) (T — t) and variance of 
a'^{T — t) as is expected for the Black-Scholes case with constant volatility [2]. 

In general for a more complicated (nonlinear) Hamiltonian such as the one given in eqn. 
(12b) for stochastic volatihty it is not possible to exactly diagonalize H and then be able 
to exactly evaluate the matrix elements of e~'^^. The Feynman path integral is an efficient 
theoretical tool for analysing such nonlinear Hamiltonians, and this is the reason for using 
the path integral formalism. 



Appendix 2 The Lagrangian 

The lagrangian L is central to the path integral formulation of quantum mechanics. We first 
find Lbs ioi the Black-Scholes case of constant volatility before tackling the more complicated 
case. For infinitesimal time e it is given by Feynman's formula 

< x\e-'^^'''\x' >= iV(e)e'^^^ 

where N{e) is a normalization constant. Since the formula (A4) is exact, we have, for 
5x — X — x' 

and with 
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We now analyze the case of stochastic voltihty. Prom eqn. (12b), the Hamiltonian is 
given by 



H 



dxdy 2dy^^2^ ^' dy 



(A6a) 



Hence we obtain using (A3) 
<x,y\(i \x ,y > = / —<x,y\e \px,Py >< Px,Py\x ,y > 

_ f°° dpx dpY ^ipx(x-x')^ipY{y-y')^-eH(x,y,px,PY) [A6b) 

and from (A4) the matrix elements of the Hamiltonian is given by 

e^ 1 ^^1 
H{x, y,Px,PY) = + (2^^ ~ ^''^^^ ^ iP^^''^PxPY + i^Py + (2^^ ~ ^^)^PY {Al) 

To perform the Gaussian integration over px and py we need the inverse and the deter- 
minant of the matrix 



Note 



and 



Let 



and 



detM = eV(l-p2) 



e(l-p2) 



-epe-^/2 1 



A — X — x' -\- er — -e^ 



B = y -y' + -i 



{All) 



{A12) 



From (A6) and (A7), on performing the Gaussian integrations, we have the Feynman relation 

1 



< a;, y\e ^^\x\ y' > 



27reVdetM 



(^13) 



where from (A6), (A8) and (AlO) 



2^^(rr7) + ^.B^ - 2f e-/^A5) + 0(e) 



{AU) 
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and finally from (All) and (A12), for 5x — x — x' and Sy — y — y', we have the negative 
definite lagrangian given by 

e"'^ fSx i Y P vn/^y , 1 



Appendix 3: The / DX - Path Integration 

The path integration over the a;(t)-variables in eqn. (19) can be done exactly, and the 
derivation is given below [9]. 
Let 

j DXq^^ (A16) 
where from (22b) the x-dependent term of the Lagrangian is given by 

Lx{i) = -^7f^{— + r - Uy, - f e^'/^(^ + - \e)Y + 0{e) {Am) 
2(1 - p^j e 2 ^ e 2 

Let 

Q = r-le^-p/^(^ + /.-ia {A17a) 

Then 

1 ^ 

^ "^TTT^^ ^ '''''' ^""^ ~ ""'-^ + 
P ) i=i 

with boundary values given by 

xq = x', xn = X {AVJc) 

We make the change of variables 

i 

Xi = Zi — Cj, dxi = dzi, i = 1,2, . . . , N — 1 {A18a) 
with boundary values 

zo = x' {Al8b) 

N 

Zn — X + e'^Cj {A18c) 
i=i 

We hence have from (A17b) and (Al8a) 

1 ^ 
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and, from (25c) and (A18a) 

Q = , n / y (Am) 

All the integrations can be performed exactly; one starts from the boundary by first 
integrating over zi , and then over Z2, ....and finally over zn-i- Consider the zi -integration. 
We have 

/+0O /7y, p~^i/2 1 
exp{--J-— [e-^^(., - z^r + e-^'(-i - -o)^]} 
-~ ^27re(l-p2) 2e(l-p^) 

e^2/2 1 1 

^M-w-r. — :5V .V. , .vJ ^2 - ^o)'} (^20) 



Ve^i + e^2 2e(l - p2) + e^2 

Repeating this procedure {N — 1) times yields 



Q = , ^= (^21a) 

/27re(l-p2)^^^^e^' 



where 



2 



N 



= ^ — M — 7f{x - x' + ey CiY (A21b) 

For the case of constant volatility = = p and e^* = o"^ = constant for all i; eqn. 
(A21b) then reduces to the Black- Scholes given in eqn. (A3). 

Appendix 4 : Generating Function Z{j,y,p) 

We evaluate 

Z{j: y,p) = j DY dtj{t)y{t)](^''^'y'e''^; {A22a) 

from (27b) we have 



with boundary condition 



y{r) = y. {A22c) 
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We first evaluate Z{j, y, y') given by the path integral in (A22a) but with the boundary 
condition 

y(0)=y',y(T)=y (A23) 

and from eqn. (A22a) 

/oo 
dy'Z{j,y,y'yy/\ (A24) 
-oo 

Define new path integration variables z{t) by 

z{t)=y{t)-y' + -{y'-y). (A25) 

T 

Note the new variables z{t) due to (A23) have boundary conditions 



Hence 



Z{j,y,y') = e'^''jDZe'- {A27) 



with 



Wo = -^(i/' -y-f^T + ^er'r + y' [ dtj{t) 



2re 

y'-y r 



T JO 



r dttjit) iA2%) 
Jo 



and 



To perform the path integral over z{t), note that from boundary condition (A26) we have 
the Fourer sine expansion 

oo 

z{t) = Yl sm{Trnt/T)zn (ASO) 

n=l 

From (A29) and (A30) we have 

2 oo oo „T- 

Sz = -—-^Y.^'^l + T.[ dtj{t)sm{7rnt/T)]zn {A31) 

n=l n=l 

The path integration over z{t) factorizes into infinitely many Gaussian integrations over 
the Zn 's and we obtain 

J DZe^^ = C J°° dzidz2dz3 ■ ■ ■ dz^oe^^ {A32) 
= C(r)e^ {A33) 
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where 

r dtdt'j(t)D(t,t')j(t') (A34) 

71^ Jo 



- r dtj dt'j{t){T-t)t'j{t') {A35) 

T Jo Jo 



Since 



D{t, t') = J2— sm(7rnt/r) sin(7rnt7r) (A36) 



n=l ^ 



'2 tt\ 



^[t'e{t - t') + te{t' -t)-—]. {A37) 
2t t 



where the step function is defined by 





1 


t > 


e{t) = < 


1/2 


t = 







t < 



{A37a) 



The normahzation function C(t) can be evaluated by first considering the discrete and 
finite version of the DZ - path integral along the lines discussed in Section III. The result is 
given [Ref 9] by 



Collecting our results we have 



qWo+W 



and performing the y' Gaussian integration in (A24) finally yields 

Z{j,y,p)^e'' {A40) 

with 

F^y r dtj{t) + {^x- \e) f dt{r - t)j{t) 
Jo I Jo 

+ e r dtj{t){T - t) f dt'j{t') + F' (AAl) 
Jo Jo 

where 

F' = IvP + [ dt{T - t)m + ^pr(/. - + r. (A42) 
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